No-go theorem for photon bunching suppression via all-unitary linear optics 
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We prove that linear optical elements such as phase shifters and beam splitters alone, cannot be used to 
simultaneously generate arbitrary single-qubit gates on a pair of qubits realized as single-photon optical modes, 
and at the same time prevent coupling to the subspace of states of two photons per mode. This constitutes a 
no-go theorem for photon bunching suppression via all-unitary linear optics. 
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Introduction. — Optical implementations of qubits have 
played an important role in quantum information science. 
Photons have an intrinsic lack of decoherence and are sim- 
ple to control by standard off-the-shelf components. Further- 
more, photonic qubits for quantum computation are particu- 
larly attractive because they can be used to interface to var- 
ious quantum communication applications jll]. Due to the 
extremely small photon-photon coupling available in existing 
materials, it was at one point believed that optical qubits could 
not be used for scalable quantum computation. However, it 
is now understood that the process of photon detection itself 
can lead to effective photon-photon nonlinearities. In partic- 
ular Knill, Laflamme, and Milburne (KLM) 12(1 launched the 
field of linear optics quantum computing (LOQC) by showing 
that deterministic single-photon sources and high-efficiency 
single-photon detectors allow the realization of scalable, prob- 
abilistic quantum computation purely with linear optical ele- 
ments. This holds in spite of the fact that using linear op- 
tics alone, the success probability of the nonlinear sign shift 
gate used in the KLM scheme cannot be improved above 1/4 
i^\M- Since the original KLM proposal, a number of authors 
have suggested various simplifications, modifications, and op- 
timizations, e.g., isl^Sl]. Although these results formally show 
that scalable quantum computing is possible, the realization is 
very demanding in practice due to the large resource overhead 
arising from the required non-deterministic photon detection 
events. This becomes particularly apparent when considering 
beamsplitter-based two-photon gates, which areprobabilistic 
in nature. In fact all optical two-qubit gates j|9t], including 
the promising non-destructive CNOT gate jldllin . fail in the 
case where more than one photon is emitted into the same op- 
tical output mode, due to photon bunching. One proposal to 
suppress such gate failure events in the polarization encoding 
uses the quantum Zeno effect, by coupling to the output light 
fields to atomic transitions il2l Il3fl . 



Here we address the natural question of whether linear op- 
tical elements alone can be used to establish quantum interfer- 
ence such that the photon-bunching effect, and thus the prob- 



ability of gate failure, can be suppressed. We show that this 
cannot be achieved in either the single-rail or dual-rail encod- 
ing paradigms, in the sense that one cannot simultaneously 
implement a linear optical two-photon entangling gate and de- 
couple the double-occupancy states 1 14]. We thus prove a no- 
go theorem for p hoton bunching suppression via all-unitary 
linear optics jlSj]. 

Linear optical elements. — In the standard circuit model of 
quantum computing any unitary transformation on n qubits 
can be decomposed as a product of gates, each of which acts 
nontrivially on at most two qubits, and is the identity on the 
other qubits jll 6]. Likewise, in the linear optics model, any 
unitary transformation on m modes can be decomposed into a 
product of "linear optical elements," each of which acts non- 
trivially on at most two modes, and is the identity on the other 
TO — 2 modes JITII . As linear optical operations we consider 
just phase shifters and beam splitters as they generate all lin- 
ear optical elements [i.e., all 2 x 2 unitaries, the group U(2)] 
when acting on the same pair of modes. Moreover, when they 
are allowed to act on overlapping pairs of modes, the linear 
optical elements generate all ni x m unitaries [i.e., the group 
U(to)] M- 

Single vs. dual-rail encoding. — Single photon optical 
qubits come primarily in two varieties; single-rail 11^ uOfl . 
where a qubit is represented by the absence or presence of a 
single photon of fixed polarization in one optical mode, and 
dual -rail 11611 . where each qubit is encoded into the presence 
of a single photon in one or the other of two spatial optical 
modes. Polarization qubits 12111 — where each qubit is repre- 
sented by two orthogonal polarization modes — are formally 
equivalent to dual-rail qubits ||9t], so we will use both inter- 
changeably. Most LOQC proposals, including the original 
KLM scheme, use dual-rail qubits, but there is considerable 
interest in the single-rail encoding scheme as well, as evi- 
denced by the numerous experiments devoted to preparation 
of arbitrary single-rail qubits, e.g., 122-26]. In the dual-rail 
encoding single-qubit operations are straightforward while 
non-trivial two-qubits are challenging and are implemented 



probabilistically y]]. The reverse is true for single -rail en- 
coding: non-trivial two qubit gates are straightforward, while 
single-qubits can be implemented probabilistically 12711 . In 
both the single and dual-rail encodings two-qubit gate failure 
due to photon bunching is a major challenge standing in the 
way of scalability. 

The subspace of interest in the single-rail encoding is 
spanned by the 6-dimensional Fock state basis AU B where 



A={l,alalalal}\0), 
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Here |0) is the global vacuum, aj, at, and rii are the creation, 
annihilation, and number operator for mode i, respectively, 
satisfying the bosonic commutation relations [a^ 
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[ai,ni\ = Ui, and [al,ni] = — a|. The subspace A is the 
computational basis of single-photon states in each of the two 
modes, while B is the subspace of two photon states in either 
one of the two modes (double-occupation states). The latter 
are the problematic states causing computational errors in the 
single-rail encoding. To define a basis for the polarization en- 
coding consider the creation operators a| , where i <E {1, 2} is 
the spatial mode index, p g {H, V} is the polarization mode 
index, and aU a] commute for all values of i and p. The ith 
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qubit basis states are then |0)j = a-^|0) and |1)- ~ o^lvl^) 
so that the computational and double-occupation (error) sub- 
spaces are, respectively. 
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S' = {al^/V2,al^/V2,a|^al^}|0), *e{l,2} (4) 

It is well known that linear optics is insufficient for gener- 
ating non-trivial (entangling) two-qubit gates in the dual-rail 
encoding. Indeed, this, together with the inherent weakness of 
non-linear photonic interactions such as the Kerr effect, is the 
raison d'etre of the KLM scheme. The reason is a straightfor- 
ward group-theoretic counting argument. Consider m — 2K 
optical modes. If dual-rail linear optics could generate an en- 
tangling gate between one pair of qubits, which could then 
be swapped to all other pairs of qubits using beam splitters, 
the set of Unear optical elements would generate the group 
U(2^). However, this contradicts the fact that the set of lin- 
ear optical elements can only generate \J{2K) 1 18112811 . This 
means that it is impossible to generate an entangling gate in 
the subspace A', irrespective of what happens in B'. However, 
this argument cannot be applied in the single-rail case, so we 
devote the remainder of this work to the no-go theorem for the 
single-rail case. 

No-go theorem for two modes. — The Hermitian beam- 
splitter generator is 



X 



- \a\ai + a|a2 j 



(5) 



and the Hermitian phase shifter generators are, for each mode, 

ni = ajoi, n2 — a\a2- (6) 



Let us also define 

r E 



( fljOi — a\a2 ) 



-(rii-ris). (7) 



Then it is easily checked using the identities in 112911 that the 
set of Hermitian operators {X, Y, Z} is closed under com- 
mutation and moreover satisfies su(2) commutation relations. 
Therefore the beam-splitter generator X and the phase shifter 
generators (combined as Z) generate U(2). The most general 
unitary evolution operator we can then construct according to 
the Euler decomposition is the following composite gate: 

C/(a,/3,7,,5,e) = e'^°'"'+'^''^^e''^ e'^'''"'+'^"'\ (8) 

where we have allowed a more general linear combination of 
the two phase shifter generators than Z, since they can be in- 
dependently tuned. 

We are interested in the (reducible) representation of the 
unitary operator ([8]) in the 6-dimensional Fock state basis 
AU B. We shall now show that if one demands that U does 
not couple A and B, then without further assumptions it is 
not possible to implement an entangling gate between the two 
single-rail qubits. In the AU B basis we have 
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C/(«,/3,7,5,e) = 
diag(l,e'",e'^e'("+«,e^'",e2'^) 

1 
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B2 



where diag denotes a diagonal matrix and where 
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(10) 



Note that {/(a, /?, 7, (5, e) consists of three blocks (1, A2, and 
B2) of dimensions 1 (vacuum), 2 (one photon), and 3 (two 
photons), respectively. Therefore repeated applications of 
the composite gate U with different values for the angles 
a, /3, 7, (5, e will still have the same triple block structure. 



Hence our construction is general in the sense that we need 
only consider a single composite gate U. 

Now, the problem is that the two-photon block includes B. 
We would like to decouple these two states from the other 
four. Clearly, this can be done by setting the matrix elements 
f/4,5, C/4,6, C/5,4 and Ue,4 (i.e., the (1, 2), (1, 3), (2, 1), (3, 1) 
elements of B2) to zero. As is clear from Eq. (fTOl l. a necessary 
and sufficient condition for this is 



sin2e = 0. 



(11) 



This is achieved whenever e is an integer multiple of tt/2. If 
we pick an even multiple 2n of tt/2 we have 

C/(a,/3,7,^,n7r) = diag(l,(-l)"e'("+'^),(-l)"e^('^+*', 



J{a+l3+-y+S) p2i(a+7) p2i{P+S) 



If we pick an odd multiple 27T, + 1 of tt/2 we have 

1 



[/(a,/3,7,<5,(2n + l)-) 



A, 



S, 



(12) 



(13) 



where 



A, 



B.s 
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_e2*("+'') 



_pi(a+l3+'r+S) 







=2J(/3+7) 







In both cases U has a block-diagonal structure where the two- 
photons per mode states are decoupled from the other four 
states. Considering just the 4-dimensional block acting on 
span(^), we have for the even case 

diag(l, (-l)"e'("+''), (_l)»e*(/3+'5), e'("+'3+T+'5)) = 
diag(l, (-l)"e'('^+'5)) ® diag(l, (-l)"e*("+T)) (14) 

i.e., a tensor product of two single-qubit phase gates. Simi- 
larly, for the odd case we have 



10 

(-l)"ie*("+'') 

(-l)"ie'('3+^) 

_g«(a+/3+7+'5) 

= SWAPdiag(l, (-l)"ze'("+^)) ® diag(l, (-l)"ie''('3+'^)), 

(15) 



where SWAP 



1000 
0010 
0100 
0001 



, a non-entangling gate. Unfortu- 



nately this means that, having decoupled span(S), the action 
of U on the A subspace is equivalent to a tensor product of 
two single-qubit phase gates, and therefore is no longer en- 
tangling. 

The reason for this no-go result is straightforward: we are 
forced to set e equal to a multiple of 7r/2 in order to prevent 



coupling to the span(,B) subspace. Once we do this we are 
left with only four independent parameters {a, (3, 7, 5), which 
are the parameters of the two phase shifters [see Eq. (|8]l]. The 
latter can only generate phase gates. To circumvent this re- 
sult would require the use of some type of nonlinear optical 
element. 

No-go theorem for M > 2 modes. — One might wonder 
whether the inclusion of additional modes allows to circum- 
vent the no-go theorem. We next show that this is not the case. 
Consider a linear optical system with 1\I modes containing 
a total of two photons, of which only the first two modes- 
supporting t wo q ubits-are the target of linear optical quantum 
computation pC!]. We shall therefore refer to these two modes 
as the "computational modes." With up to two photons per 
mode these comprise the computational basis A ^ and the 
undesired double-occupation states B (|2]i. The other 1\I — 2 
ancillary modes are used to try to remove the span(i3) sub- 
space of the computational modes. The relevant basis states 
are therefore AUB for modes j = 1 and 2, and the set of Fock 
states 



(16) 



for modes j ~ 3, . . . , Af, where each occupation number 
Uj e {0, 1, 2} v.? e {1, . . . , M}, and J^^L^ "i < 2. The 
complete 2-photon Hilbert space is then H = span[(^ US)® 
C]. 

Now let V be some arbitrary linear unitary optical transfor- 
mation on all, or some subset, of the M modes. For example, 
V could be a product of phase shifters and beam splitters cou- 
pling any pair of modes. By linearity, V has the following 
matrix representation; 




VaW"^ 



M 

E 



v,,a 



e{i,...,M}. 



(17) 



Note that this matrix representation is in the basis of bosonic 
creation operators, which differs from the Fock state basis 
spanning H. We need both of these representations for our 
proof. Whatever V is, we require that after it is applied the 
computational Fock state basis states A are decoupled from 
the rest. In other words, while V may be constructed from lin- 
ear optical elements which couple different states, including 
the computational basis and all other states, at the end there 
may be no transitions between the computational subspace 
span(^) to the rest of the Fock space H' = H \ span(^). 
Formally, this means that we require that for all computational 
basis states \tfja) G A and for all other basis states {ipp) G H', 



{ijo.\V\^^)^0. 



(18) 



This requirement states that V has a block-diagonal struc- 
ture in the Fock state basis; V ~ Va ffi V-h', where Va is 
a 4 X 4 matrix that is supposed to implement an entangling 
operation inside the computational basis, and Vw is some 



dimH' X dim'H' linear optical matrix. We shall now show 
that by combining this Fock space requirement with the lin- 
earity of V (as expressed in the basis of creation operators), 
the problem reduces to the one we already studied in the two- 
mode case, namely, Vh' = Vg ® Vw, where V13 is a 2 x 2 
matrix acting exclusively on the B states, and Vh" is some 
linear optical matrix acting in the H' \ span(i3) subspace. 
Because of this block-diagonal structure of V, in particular 
V4 © Vis, the conclusion from the Af = 2 analysis, that V_a is 
non-entangling, will apply directly, i.e., V^ will have to be in 
the form of either Eq. (fl4b or Eq. ( fTSl l. 

To prove this, consider the matrix element of V between 
a single -photon computational basis state a[\0) and a single- 
photon state in some mode j £ {3, . . . , Af }, which must van- 
ish by our no-leakage requirement ( fTSb : 



vanishes: 



= (0|ajya||0) = {0\ajVa\v''V\0) 

M 

= J2vik{0\ajalv\0) 



(19a) 
(19b) 



k=l 
2 



M 



J2vik{0\aia,V\0) +J2^ik{0\[ala, + 6,k)v\0 

A;=l fe=3 

vi,{0\V\0); je{3,...,A/}, 



(19c) 
(19d) 



where to arrive at ( |19bt we used the linearity condition 
Eq. iT7\ . to arrive at ( |19c| ) we used the commutation rela- 
tion [flj, a|,] = Sjk, and to arrive at the last line we used 
(OjaJ. = 0. Now, because y is a product of linear optical 
elements, and each such element can be written as cxp(— iil), 
where fl is linear combination of X, Y and Z [Eqs. (|5]l and 
(|7]i, with the appropriate mode index adjustment], we have 
r^lO) = 0, so that V\0) = |0), and hence (0|y|0) == 1. We 
can thus conclude that vij = 0. A similar calculation leads to 
Vji = V2j = Vj2 = 0. These constraints state that, subject to 
Eq. ( fTSb . V does not couple modes 1, 2 and the ancilla modes, 
i.e., it has a block-diagonal structure in the creation operator 
basis. 

Next, let us show that this result means that V cannot cou- 
ple the double-occupation states B to the ancilla modes, i.e., 
that Vh' = Vb®V-u"- Note first that because \^ is a linear op- 
tical unitary, it commutes with the number operator, i.e., it pre- 
serves the mode occupation number, and hence its only non- 
vanishing matrix elements are between Fock states with equal 
occupation numbers. For this reason we need only consider 
matrix elements of V between elements of B and two-photon 
states in C, e.g., {G\a\V a\a}-\Q) , where i,j e {3,..., A/} 
(i = j OT i y^ j). Let us show that this matrix element also 



{0\alValal\0) = (0\alValvWalvW\0) (20a) 

M 

= J2 y:kVji{0\alalalv\0) (20b) 



k,l = l 



M 



= ViiVji{0\ala\a\\0) +'^VikVji{0\alala\\0) 

k=2 
M M 

1=2 fe,/=2 

(20c) 
= 0; i,je{3,...,A/} (20d) 



where to arrive at ( |20bt we used Eq. ( fTTI i again, to arrive at 
( I20cb we used V\0) = |0) along with the bosonic commuta- 
tion relation [a^, a A = 5ij, and to arrive at ( I20dl ) we used the 
previously derived vij = Wji = V2j = Vj2 = for j > 3 for 
the first three terms in ( I20cl ) (the second and third terms also 
vanish due to (0|aj, = {0\aj = 0), and ai |0) = for the last 
term. A similar calculation shows that (0|a2V^a]^aJ||0) = 
for i,j& {3, . . . , A/}. Thus, we have shown that in the Fock 
state basis, V has the claimed block-diagonal matrix repre- 
sentation, namely V ^ Va®Vb ® Vw ■ The M = 2 result 
thus applies, and V is equivalent to a product of two single- 
qubit phase gates inside the computational subspace span(^). 
We conclude that there is no advantage to considering a larger 
number of modes. This result is a consequence of the assump- 
tion (T% that y is a unitary linear optical transformation, and 
the requirement ( fTSI l that the computational basis states are 
decoupled from the rest of the Fock space. 

Conclusions. — We have shown that using linear optics 
alone it is not possible to suppress the photon bunching effect. 
This means that the probability of a successful gate operation 
using passive optical elements, such as non-polarizing and po- 
larizing beamsplitters, cannot be improved. This result agrees 
with previous theoretical work IJ, ISlf showing that passive 
linear optics does not involve particle interactions other than 
those imposed by statistics, and can be understood in terms of 
classical wave mechanics. Therefore additional nonUnear op- 
erations such as photon-detection or absorption 11211 . the Kerr 
effect, or light squeezing, are needed in order to enable uni- 
versal linear optical quantum computing. 
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